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ABSTRACT 

We present 1" resolution CARMA observations of the 3 mm continuum and 95 GHz CH3OH masers 
toward 14 candidate high mass protostellar objects (HMPOs). Dust continuum emission is detected 
toward seven HMPOs, and CH3OH masers toward 5 sources. The 3 mm continuum sources have 
diameters < 2 x 10^ AU, masses between 21 and 1200 M0, and volume densities > lO^cm"'^. Most of 
the 3 mm continuum sources are spatially coincident with compact H II regions and/or water masers, 
and are presumed to be formation sites of massive stars. A strong correlation exists between the 
presence of 3mm continuum emission, 22 GHz H2O masers, and 95 GHz CH3OH masers. However, 
no 3 mm continuum emission is detected toward ultracompact H II regions lacking maser emission. 
These results are consistent with the hypothesis that 22 GHz H2O masers and CH3OH masers are 
signposts of an early phase in the evolution of an HMPO before an expanding H II region destroys 
the accretion disk. 

Subject headings: stars; formation - masers 



1. INTRODUCTION 

The evolutionary sequence of low mass star formation 
is broadly understood. Gravity in a prestellar core even- 
tually overcomes the thermal, turbulent, magnetic and 
rotational support to form a protostar. This protostar 
is surrounded by an envelope and accretion disk, and 
launches a bipolar outflow. As the protostar accretes 
mass, the surrounding envelope and disk are cleared by 
the winds and outflows until no more m ass is left to ac- 
crete (see, e.g., iMcKee fc Ostrikeri[2007l) . 

High mass protostars share many of the same observa- 
tional characteristics as th eir lo w mass counterparts: out- 
flows (|Zhang et al.' '200 11) . iets |Caratti O Garatti et all 
I2OO8D . envelope s (Mohn aji et al.lll996f ) and possibly disks 
(jCesaroni et al.ll2005il) . However, the elevated radiation 
environment that accompanies high mass stars raises the 
question as to what degree high mass star formation is 
merely a scaled version of low mass star formation. An- 
swering that question is not straightforward, partly be- 
cause high mass star forming regions are less common 
than the low mass counterparts, and tend to be more dis- 
tant. The greater distance, coupled with the propensity 
of HMPOs to form in clusters and the disruptive influ- 
ence of jets, outflows and H II regions, complicates inter - 
pretation of the observations (jZinnecker fc Yor"kil2007l ). 
Therefore, to study the initial conditions of high mass 
star formation, it is important to identify young HMPOs 
before their surroundings have been altered by energetic 
processes. 

Surveys for high mass protostars have been con- 
ducted over the years at a variety of wavelengths. 
Optical and radio continuum surveys of the galac- 
tic plane have identifled H II regions in which mas- 
sive star format i on is already well underway (e.g. 
Marsalko ya) 119741: 'Ahenho ff et all Il979t iHavnes et all 
1979). W ood fc Ch urchwcU ([l989l ) used the IRAS cat- 
alog to search for far-infrared sources associated with 
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ultracompact H II (UCHII) regions and high mass pro- 
tostars embedded in dense cores. Massive star forma- 
tion is thought to be associated with Infrared Dark 
Clouds (IRDCs), which are seen in absorption against 
the diffuse mid-infrared background of t he galactic plane 
(|Perault et al.l 119961 : lEgan et al.l [TUl). The masses 
of IRDCs suggest that they may form stellar clusters 
(|Rathborne et al.l [2006. 2007), and IRDCs are known 
to ho st both protostars and starless cores (jWang et al.l 
l2006l ). More recently, surveys with the Spitzer Space 
Telescope such as GLIMPSE (jBeniamin et all l2003f) 
have id entifled HMPOs based on mid-infrare d colors 
(jKumar & Grave 2007; Cvgano wski et al.ll2008l ). 

High mass star formation also coincides with a va- 
riety of masers, as demonstrated by targeted sur- 
veys of known HM POs and blind surveys o f the 
galactic p lane (e.g. Szvm czak. Hrvnek. fc Kui l2000t 
iPandian et al. 2007). Intcrferometric observations have 
shown that water masers (at 22.2 GHz) are associ- 
ated with J-shocks fr om jets and winds ejected by 
massive YSOs ( Torrcllcset_aD 120031: iGoddi et al.l l2005t 
iMoscadelh et al.ii2005i) . Methanol masers (at 6.7 GHz) 
are less understood, but are thought to trace the 
shock where a j et meets the ambient circumstellar gas 
(|De Buizeij|2003[ ) an d/or the innermos t portions of a disk 
around a protostar (|Edris et al.ll2005f ). Hydroxyl maser 
emission (at 1.7 GHz) may also come from a disk around 
an accreting protostar (Edris et al. 2005) , an d is also as- 
socia ted with both H II and UCHII regions (|Edris et al.l 
I2OOI . 

Although there are many tracers of high mass star for- 
mation (masers, radio continuum, infrared objects), they 
are not always spatially coincident or even found within 
the same object. Over the years, various observations 
of HMPOs have been unified u nder a common evolu- 
tionary sequence (see reviews by Garav fc Lizanol 
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'Menten. P illai. fc Wvrowskil l2005t iBeuther et al 
Zinncckcr fc Yorkd I2007D . IRDCs are thought to host 
the youngest protostars, perhaps surrounded by accre- 
tion disks and driving jets and outflows. The disk, jet, 
and outflow system are conducive for the formation of 
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H2O, CH3OH, and possibly OH maser emission. This 
evolutionary stage is also distinguished by dust emission 
from the disk and envelope around the protostar. The de- 
velopment of an H II region is suppressed initi ally by the 
continued accretion of material onto the star (W almslevI 
flQQSi ). As the protostar increases in mass, temperature, 
and luminosity, an H II region forms that is initially grav- 
itation ally bound to the star and remains compact in size 
(|Ketd 2002). Eventually the protostar becomes massive 
enough to produce an H II regio n tha t escapes the stel- 
lar gravitational potential (Ketoi l2007l ). The H II region 
expands to engulf the accretion disk and envelope, which 
destroys the source of the masers, jets, and outflows 
(iLo etal.1 119751: IGenzel fc Downed IiQTTI: iCodella fcFbllil 
Il995t ICodella et al.l |2004[ ). This evolutionarv sequence 
has provided a framework to explain a multitude of obser- 
vations (e.g iMolinari et al. 1996: Beuther et al.ll2002U 
iMinier et al.ll2005Hvan der Tak fc Mentenll2005f) . 

An observational prediction of this evolutionary se- 
quence is that compact dust continuum emission from an 
accretion disk or envelo pe will be present in H MPOs con- 
taining maser emission (ICodella et al.l[200l . At a more 
advanced stage when the H II region has grown, maser 
emission will be suppressed and dust emission from the 
disk and envelope should be at a much reduced level. Sin- 
gle dish maps of the dust emis sion exist for a number of 
HMPOs (jBeuther et al.ll2002af) . but no clear correlation 
is present between the dust continuum luminosity and 
the presence of maser emission. However, the resolution 
of these observations (11") is insufficient to distinguish 
compact emission around HMPOs from the ambient 
dense core. Interferometric observations have achieved 
significantly higher resolution (~5"), but they cover a 
limited range of presumed evolutionary states and are 
insufficient to test the proposed evolutionary sequence 
^ohnari ct al. 2002; Beltran et al. 2004; Cesaroni et al. 
[20051: iZhang et al.ll2007l : IReid fc Matthewdl2008D . 

In this paper, we present high resolution (1") CARMA 
observations of the dust continuum and 95 GHz CH3OH 
maser emission toward 14 HMPOs to test the evolution- 
ary sequence for massive stars described above. The 
HMPO sample was selected to cover the primary evo- 
lutionary stages in this model. With the order of mag- 
nitude improvement in angular resolution over previous 
single dish surveys, we can isolate any compact dust 
emission originating from a disk or envelope associated 
with an HMPO. In iJH we describe the sample of HMPOs 
in our survey, and the available ancillary information 
regarding their evolutionary state. The new CARMA 
observations are described in ^3] and the observational 
results are presented in SJ?) In ^Sl we compare our re- 
sults with the proposed evolutionary sequence for mas- 
sive stars. 

2. SAMPLE 

We drew a sample of 14 high-mass protostella r can- 
didates from t he survey _ofJSridharaneLalJ (12002). The 
HMPOs in the iSridharan et al,, (.200211 list have I RAS co\- 
ors of UCHII regions (I Wood fc Churchwel]|[l989l) are de - 
tected in the CS J=2-l survev of lBronfman et al.l(|1996D . 
are bright at FIR wavelengths (Fgo > 90 Jy, Fioo > 500 
Jy), are north of —20° declination, and are not detected 
ab ove 25 mJy in the 5 GHz radio continuum su rveys 
of iGregorv fc CondonI (fl99l[l . iWright et all ([1991 and 



iGriffith et all (fT99l . 

The properties of the HMPOs selected for this study 
are summarized in Table[TJ The HMPOs have a bolomet- 
ric luminosity between 10^'^ and 10^ ° L©. The distances 
to the sources range from 1.7 to 10.5 kpc, where the 
near-far distance ambiguity has been r esolved through 
a vari ety of methods as described in ISridharan et al.l 
()2002D . The dense cores associated with the HMPOs 
have masses of ~4OO-1OOOOM0 contained within a ~1 pc 
region. We adopted the velocity and association with cm- 
continuum emission, H2 O masers and CH3 OH masers 
for each source given in ISridharan et~al] (|200l . The 

3.6 cm radio continuum emission was observed with the 
VLA in B-array with a sensitivity of cr ^1 mJy, and 
the coordinates of the detections were kindly provided 
by H. Beuther (private communication). The 1.2mm 
continuum fluxes and positions were measured with 
MAMB O on the IR AM 30m telescope at a resolution 
of 11" ([Beuther et a l. 2002a i). The assoc iation with OH 
masers is presented in Edri s et al.l ([20071 ). and were de- 
tected with the Nancay radio telescope and the GBT. 

The 14 HMPOs that we observed were divided into 
three categories depending on the presence or absence 
of cm-continuum emission indicative of an UCHH re- 
gion, 6.7 GHz CH3OH masers, or 22 GHz H2O masers. 
In order of increasing age according to the proposed 
evolutionarv sequence ([Lo et al.lll975l : iGenzel fc Downed 
[19771: ICodella fc Fellilll995D . our sample contains (1) five 
sources with (CH3OH or H2O) maser emission but no 
cm-continuum, (2) five sources with both cm-continuum 
and maser emission, and (3) four sources with cm- 
continuum but no maser emission. For each of these 
three categories, the range of distances (2.9 — 9.5 kpc, 

1.7 — 10.5 kpc, 2.6 — 5.5 kpc, respectively) and luminosi- 
ties (lO" ! - IO^ OLq, 103-9 - 104-9 Lq, lO"-! - IO^-^Lq) 
sampled are similar. We centered each CARMA map 
on the position of the peak 1.2 mm continuum emis- 
sion d erived from 11" resolution single- dish observations 
(Beut her et all l2002al) . Note that in ISridharan et all 
(|2002f l the source IRAS 18345-0641 is associated with 
both maser and cm-continuum emission. For the pur- 
poses of this survey, we consider IRAS 18345—0641 to 
be a maser-only HMPO because the masers and 1.2 mm 
continuum peak are coincident to within 1", but the cm- 
continuum emission is offset from the maser position by 
44", or 2 pc. 

3. CARMA OBSERVATIONS 

Continuum and spectral line observations in the 3 mm 
window were obtained with CARMA (Combined Array 
for Research in Millimeter- wave Astronomy). CARMA 
is a 15 element interferometer consisting of nine 6.1 me- 
ter antennas and six 10.4 meter antennas. Data were 
taken in the CARMA B-array configuration between 
2007 December 24 and 2008 February 11. The pro- 
jected baselines in this configuration range from 65m 
to 800m. The CARMA correlator records signals in 
three separate bands, each with an upper and lower side- 
band. Two bands were configured for maximum band- 
width (468 MHz with 15 channels per band) to observe 
continuum emission, providing a total continuum band- 
width of 1.87 GHz. One band was conflgured with 31 
MHz bandwidth across 63 channels (with a resolution of 
0.488 MHz or 1.6 kms^^ per channel) to observe the 
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hyperfine lines of CH3CN (5-4) in the lower sideband 
(from 92.959 to 92.987 GHz) and CH3OH (8o,8-7i.7) in 
the upper sideband (at 95.2 GHz). The FWHM of the 
synthesized beam with natural weighting of the visibil- 
ities is ~1" (see Table [2|), and the largest angular scale 
that can be accurately imaged in the maps is ■~5". The 
half-power beam width of the 10.4m antennas is 73" at 
the observed frequencies. 

The observing sequence was to integrate on a primary 
phase calibrator for 2.5 minutes, an HMPO for 7 min- 
utes, and a secondary phase calibrator for 2.5 minutes. 
This cycle was repeated for ~3 hours for each HMPO. 
Calibration and imagin g were done using the MIRIAD 
data reduction package (jSault et"an [l995). The sources 
observed, along with the phase center, calibrators and 
the noise in the resultant continuum maps and spectra 
are listed in Table [3 The secondary calibrator was used 
to test the accuracy of the phase transfer from the pri- 
mary calibrator. A passband calibrator (3C454.3) was 
observed for fifteen minutes in each set of observations, 
and radio pointing was performed every two hours there- 
after. Absolute flux calibration was accomplished using 
Neptune, Uranus and MWC349 as primary flux calibra- 
tors and 3C454.3 as a secondary flux calibrator. Based on 
the repeatability of the quasar fluxes, the estimated ran- 
dom uncertainty in the measured source fluxes is cr 5%. 
Because we spent only ~3 hours on each of the four- 
teen HMPOs we had sufficient signal to noise to detect 
CH3CN in only one source, IRAS 20126-^4104, and the 
CH3CN observations are not discussed further. 

4. RESULTS 

Here we present the CARMA 3 mm continuum and 
95 GHz CH3OH maps and compare them with maps of 
6.7 GHz CH3OH and 22 GHz H2O masers and radio 
continuum emission taken from the literature. We derive 
masses for the 3 mm detections, and present evidence for 
outflows eminating from these sources. 

4.1. 3 mm continuum 

We detect 3 mm continuum emission toward seven of 
the fourteen HMPOs at the 3cr level or greater. Con- 
tour plots of the 3 mm flux density for the detections are 
presented in Figure [U and plots of the non-detections 
are shown in Figure [2l The sources IRAS 23033+5951 
and IRAS 19217+1651 are resolved into three and two 
components, respectively. None of the detected sources 
are more than a factor of ~2 larger than the ~1" beam 
of our maps. The secondary phase calibrators, assumed 
to be point sources, have the same range of sizes as the 
detected HMPOs, so all of the 3 mm detections are con- 
sistent with point sources at 1" resolution blurred by 
seeing effects. The CARMA observations are sensitive 
to angular scales as large as 5", and resolve out the 
majority of the extended emission associated with the 
star-forming regions s een in the 1.2 mm single dish maps 
(jBeuther et al.ll2002aD . 

The positions and measured flux densities of the 3 mm 
continuum sources are listed in Table [H The brightest 
3 mm continuum source detected in each map is located 
within 6" of the centroid of the larger scale core traced 
by single dish 1.2 mm co ntinuum emission as measured 
by iBeiather et al.1 ()2002aD . Given the 11" resolution of 
the 1.2 mm observations and the difficulty in finding the 



exact center of non-symmetric intensity distributions, the 
offsets between the 3 mm and 1.2 mm positions are most 
likely not significant. 

4.2. Masers 

We detect 95 GHz CH3OH emission toward five HM- 
POs. Two of the sources (IRAS 23033+5951 and IRAS 
23151+5912) have a single CH3OH emission detection, 
and three (IRAS 18517+0437, IRAS 20126+4104 and 
IRAS 23139+5939) have multiple detections. Spectra for 
each of the detections are shown in Figure [3l and spec- 
tral parameters are summarized in Table [J] The 95 GHz 
CH3OH lines are bright (with brightness temperatures 
around a few hundred Kelvin), narrow (unresolved at 
1.6 kms^ resolution) and spatially unresolved. There- 
fore the inferred brightness temperatures are lower limits 
and the lines are likely to be masing. 

The five HMPOs that have 95 GHz CH3OH masers are 
also detected in the 3 mm continuum. The location of the 
presumed masers are marked on the 3 mm contour maps 
in Figure [TJ Two 95 GHz CH3OH masers are coincident 
with the 3 mm continuum emission to within 0.6" (or 
< 2000-3000 AU), while the other nine masers are sepa- 
rated by ^10"* AU from the 3 mm continuum peaks (see 
Table |4]). The projected spatial offsets are consistent 
with the expectation that the 95 GHz masers are colli- 
sionally pumped at the i nterface between outfi ows and 
the interstellar medium (|Minier fc Booth|[200l) . All of 
the CH3OH lines are found within 2 kms~^ of the sys- 
temic velocity of the HMPO, implying that the masers 
form in the colliding gas. Our finding that the maser 
line is found at nearly the same velocit y as the ambient 
mediu m is consistent with the results of lMinier fc BootU 
(|2002f) . who find that this maser line is often, but not 
always, at the systemic velocity of the HMPO. 

The location of known centimeter continuum sources, 
H2O masers, and CH3OH masers are marked on the 
contour maps in Figures [1] and O Six of the seven 
CARMA 3 mm continuum detections (all except IRAS 
18517+0437) were included in an interferometric sur- 
vey of 22GH z H2O and 6.7 GHz CH3OH masers by 
iBeuther et al] (l2002bf ). Water masers have been iden- 
tified toward all six cores detected with CARMA, and 
three of these cores also have 6.7 GHz CH3OH masers. 
For HMPOs with both 22 GHz H2O and 6.7 GHz CH3OH 
masers, the two types of masers are spatially coincident 
to within 1-2". Both types of masers also fall within 
~1" of our 3 mm continuum positions, as shown in Fig- 
ure [T] Of the multiple 3 mm components resolved in 
IRAS 19217+1651 and IRAS 23033+5951, one is coin- 
cident with H2O masers in each of the HMPOs. 

4.3. Radio Continuum 

Four of ten 3 mm continuum sources are located within 
1" of a 3.6 cm radio continuum source. The 3 mm con- 
tinuum emission for these objects may originate from 
free-free or thermal dust emission. If we assume that 
the sources in our sample are typical of other UCHII re- 
gions, the free-free emission will be optically thick (i.e. 
5*1/ oc ir^) for V < IQ — 15 GHz, and optically thin 
(Su PC t/~°-^) for higher frequencies (iHunter et al.lll998l: 
Turner fc BecS 120041 : IZapata et al] l2008t iHunter et all 
20081) . For a "tvpical" UCHII radio continuum spectrum, 
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the expected 3 mm continuum flux density from free-free 
emission will be less than a factor of 2.7 higher than the 
3.6 cm flux density. We estimate the relative contribu- 
tions of free-free and dust emission to the observed 3 mm 
flux density using available data from the literature. 

The s ource IRAS 19217+ 1651 was observed with the 
VLA bv lGarav et al.l (|2007f) . who measured a 7mm flux 
density of 50.4±0.16 mJy towards the brighter of the two 
CARMA detections at 3 mm. The 3mm flux density of 
the brighter component is 66 ± 1 mJy; an additional cal- 
ibration uncertainty of ~10% should be included in both 
the 7 mm and 3 mm flux density measurements. The 
7 mm det ection is consistent with optically thin free-free 
emission (jGarav et al.ll2007l ). and implies that the free- 
free contribution to the 3 mm flux density is ~46 mJy. 
The HMPO IRAS 20126-1-4104 w as observed w i th th e 
VLA between 20cm and 7mm bv iHofner et all (|2007D . 
who conclude that free-free emission should contribute 
~4 mJy at 3 mm, which is a small fraction of the 33 ± 2 
mJy detected with CARMA. Our CARMA observations 
detect three regions of 3 mm emission towards the HMPO 
IRAS 23033-1-5951, of which the northe rnmost is associ- 
ated with 3.6 cm emission of 1.7 mJy (jSridharan et al.l 
l2002f ). For a "typical" UCHII radio continuum spectrum, 
the free-free emission would contribute less than 4.6 mJy 
of the 8 ± 1 mJy detected with CARMA. The HMPO 
IRAS 23139-^5939 has a flux density of 1.4 mJy at 3.6 cm 
(jSridharan et all 120021 ) . so we assume that the free-free 
emission contributes no more than 3.8 mJy of the 12 ± 2 
mJy detected with CARMA at 3 mm. For regions with- 
out 3.6 cm emission, we assume that dust accounts for 
all of the observed 3 mm flux density. 

4.4. Masses 

We derive the total (gas and dust) mass for each of 
the CARMA continuum detections from the observed 
3 mm flux density, assuming optically thin emission from 
isothermal dust, using the equation 



M = 



(1) 



where is the integrated flux density, fd is the frac- 
tion of 3 mm emission from dust, D is the distance, 
is the dust emissivity and B^(T) is the Planck func- 
tion. The value of the dust emissivity is extrapolated 
from a value of Ko = 0.1 cm^ g~^ at 250 /zm assumin g a 
power law spectral index of /3 = 2 (jHildebrandl 1 1 983l i.e. 
Hi/ — Ko{v/vo)f^). For 95 GHz continuum observations, 
the mass can be derived as 



M 



^'^^ (mJy) (kpc 



T 
50 K 



(2) 



where we use the Rayleigh-Jeans approximation to the 
Planck function and assume a temperature of 50 K. Since 
the 3 mm continuum emitting regions are not resolved 
by our observations, a lower limit on the volume density 
is derived assuming spherical geometry with a diameter 
equal to the resolution and a mean molecular weight per 
free particle of 2.37. Estimated masses and densities for 
the individual 3 mm continuum sources are summarized 
in Table [31 The dominant source of uncertainty in the 
derived masses is th e dust emissivity, which is uncertain 
by a factor of a few (jHildebrandl 1 1 983f ) . 



The masses of the sources detected in the CARMA 
3 mm continuum maps range between 21 and 1200 M© 
with volume densities in excess of 10® cm^'^. These 
masses and densities can be compared to the large 
scale core traced by si ngle dish 1.2 mm continuum maps 
(jBeuther et al.|[2002aD . Assuming that the large scale 
core and the compact clumps have the same tempera- 
ture and dust emissivity, we flnd that the compact 3 mm 
sources contain '^20% of the total core mass, and have 
volume densities that are larger by a factor of ^^1000 
than the average core density. 

4.5. Outflows 

All 14 HMPO can didates in our sample were observed 
bv lSridharan et al.l (2002) with the CSO or IRAM 30 m 
telescopes to measure the CO (2-1) line profile. All were 
seen to have broad-line wings consistent with bipolar out- 
flows, suggesting that outflow activity is common in these 
high-m ass star forming regions. The ISridharan et al.l 
(j2002f ) survey also searched for SiO (2-1) emission, an- 
other common outflow tracer, from 12 of the 14 HMPO 
candidates in our sample (all except IRAS 18517-1-0437 
and IRAS 20332-1-4124), resulting in seven detections. Of 
the seven HMPO candidates with SiO (2-1) emission, six 
are associated with 3 mm continuum emission detected 
by CARMA, and none of the HMPO candidates without 
SiO (2-1) emission were detected with CARMA. Fur- 
thermore, of the five HMPO candidates that have nearby 
95 GHz CH3OH maser emission, four are also associated 
with SiO (2-1) e mission and the fifth w as not included 
in the survey of ISridharan et all (j2002l ). However, the 
resolution of the single-dish observations (10-30") is not 
sufficient to determine whether or not the HMPOs de- 
tected with CARMA are driving any these outfiows. 

Interferometric observations have reported detections 
of outflows from 5 of the 14 HMPOs in our CARMA 
su rvey: 

1) iZhang et all (j2007l ) observed SiO (2-1) at 5" resolu- 
tion towards the HMPO IRAS 18566-^0408, and found 
that the emission traces a collimated outflow, originat- 
ing from the position of the CARMA 3 mm continuum 
detection. 

2) SiO (2-1) and CO (2-1) emission fro m the HMPO 
IRAS 19217-1-1651 were observed by iBeuther et all 
(j2004( ) at 2-6" resolution. There is a bipolar outflow seen 
in both transitions, originating from the position of the 
CARMA 3 mm continuum source. 

3) IRA S 20126-H4104 has bee n observed at high resolu- 
tion bv I Shepherd et all (j200(]D . who report the presence 
of a N-S oriented molecular outflow and a NW-SE ori- 
ented jet centered on the position of the CARMA contin- 
uum source. Shocks identified by H2 and [S II] emission 
trace the jet, and are coincident with the positions of the 
95 GHz CH3OH masers detected by CARMA, in agree- 
ment with the hypothesis that these maser lines are col- 
lisionally pumped at the interface between jets/outflows 
and the ambient medium around the protostar. 

4) CARMA found three 3 mm continuum peaks towards 
IRAS 23033-1-5951, one associated with 3.6 cm emis- 
sion, one with a H2O maser, and one with 95 GHz 
C H3OH maser emiss i on. T his region has been studied 
bv iReid fc Matthewd (j2008l ). who study the outflows in 
this region at ^^6" resolution. There are at least two 
outflows in this region, one identifled by HCO+ (1^0) 
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emission, and one identified by CH3OH (2-1) and SiO 
(2-1) emission. The CARMA 3 mm continuum source 
coincident with the H2O maser is most hkely to be driv- 
ing the SiO outflow, and it is unclear what drives the 
HCO+ outflow or is exciting the 95 GHz CH3OH maser. 
5) The 3 mm continuum source detected with CARMA 
in IRAS 23151+5912 is c oincident with "Peak 1" in the 
875 /xm continuum map of lBeuther et al.l (|2007D and with 
the 3.4 ni r n and 1.3 mm continuum source detected by 
iQiu et al.l(|2007f) . There are at least two outflows, identi- 
fied by SiO (2-1) and (8-7) emission, in t his region, one 
orien ted roughly W-E and the other N-S (jBeuther et al.l 
I2007t lOiu et al.l[200l . The outflows are likely driven by 
two protostars, and it is not clear if the 3 mm continuum 
emission detected with CARMA corresponds to one or 
both protostars, or if the 95 GHz CH3OH maser is re- 
lated to the interaction between the W-E outflow and 
the ambient medium. 

5. DISCUSSION 

One proposed evolutionary sequence for HMPOs is 
that they originate from low or intermediate mass pro- 
tostars that are accreting material from t he surrounding 
core, presumably via a circumstellar disk (jBeuther et al.l 
l2007f) . As the protostar increases in mass, an H II 
region ultimately forms that eventually expands and 
halts accretion. While it is unclear if H2O or 6.7 GHz 
CH3OH masers originate from the disk or outflow, the 
observational evidence suggests these masers trace an 
early stage of evoluti on before the H II region has ex- 
panded signifl cantly (i Lo et al.) 119751: 



fl977l:lWalsh et al. 1998: lCodefla et al.l 
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20041). In this sce- 



nario, a young HMPO with an accretion disk would be 
identified by its water and/or methanol masers and mm- 
wavelength thermal emission from the dust around the 
protostar. When the HMPO is "intermediate aged" , an 
UCHII region detectable with cm-wavelength continuum 
emission would grow but will not immediately disrupt 
the accretion, so masers and dust thermal emission would 
also be present. As the UCHII region expands, it will de- 
stroy the disk and cut off the masers. We now evaluate 
whether or not the CARMA observations are consistent 
with this evolutionary sequence. 

First, we discuss if the compact sources detected in the 
CARMA 3 mm continuum maps signify the presence of a 
high mass protostar. The masses of the compact sources 
range between 21 and 1200 M© contained within a region 
< 10,000 AU in diameter. Thus these compact 3 mm con- 
tinuum sources are substantially denser and more mas- 
sive than the dense cores found in nearby clouds known 
to be forming low mass stars. Even if we assume only 
10% of the core mass forms into a star, a mid B-type 
star can potentially form in these cores. Moreover, four 
of the 3 mm continuum sources are spatially coincident 
with compact H II regions, confirming directly the pres- 
ence of a high mass star. Therefore, we assume that all 
of the compact 3 mm continuum sources represent dense 
dust and gas surrounding a high mass star or proto- 
star. We cannot determine if the m aterial is distributed 
in a disk or compac t envelope (e.g. [B eltran "e ran[200l 
ICesaroni et"alll2005t iReid fc Matthe ws 2008), however, 
as the dust continuum emission is unresolved by these 
observations. 

Since 3 mm continuum emission is detected toward 7 



of the 14 HMPOs in our survey we investigate whether 
the 50% detection rate is set by physical differences be- 
tween the high mass cores or by selection bias. A bias 
may be introduced since the HMPOs span a factor of 
'^6 in distance and a factor of ^10 in luminosity, but 
the sensitivity in each of the CARMA maps is approxi- 
mately the same. In the top panel of Figure [H we plot 
the noise in the CARMA 3 mm continuum maps as a 
function of the distance to the HMPO for sources with 
(black circles) and without (gray circles) CARMA con- 
tinuum detections. No clear differences in the noise char- 
acteristics or source distances are found between the two 
samples. The middle panel shows the core masses derived 
from single dish 1.2 mm continuum maps versus the bolo- 
metric luminosity of the core. Again, sources with and 
without 3 mm CARMA continuum detections share sim- 
ilar distributions, indicating that the CARMA observa- 
tions did not preferentially detect cores with the highest 
mass or luminosity. Finally, the bottom panel in Figure 
[4] shows how the peak flux density in the MAMBO single 
dish 1.2 mm continuum maps varies with the integrated 
1.2 mm flux density. We detect (at the > 3a level) 3 mm 
continuum emission from all seven of the HMPOs with 
1.2 mm peak fluxes above 300 mJy/(ll")^, while none of 
the seven HMPOs below that threshold are detected. 

Is there an observational bias for detecting only those 
sources with the highest peak surface brightnesses at 
1.2 mm, or is there an intrinsic difference between those 
cores above and below the 300 mJy/beam threshold? 
The peak 1.2 mm flux densities in our sample are all > 
100 mJy/beam. If the 1.2 mm flux density in the 11" 
single dish beam originates from an unresolved disk or 
envelope, we would expect to observe flux densities of 
>4 mJy/beam at 3 mm assuming an emissivity spec- 
tral index of (3 = 1, wh ich is appropriate for disks 
(Beckwith fc Sargentlfl99ll ). The noise in the 3 mm con- 
tinuum maps is approximately 1 mJy/beam, and all 14 
HMPO candidates would be detectable at > 3<j level if 
the central emission were compact. 

Given that all 14 HMPO candidates could have been 
detected with CARMA, we consider whether differences 
in the CARMA 3 mm continuum detection rate versus 
peak 1.2 mm flux density could be explained by the 
presence or absence of an unresolved disk or envelope. 
The median integrated 3 mm flux density detected with 
CARMA is lOmJy, which, assuming an emissivity spec- 
tral index of P — 1, would contribute ^180mJy/beam of 
emission at 1.2 mm. This is approximately equal to the 
difference between the peak 1.2 mm flux densities of the 
faintest HMPOs detected by CARMA (- 400 mJy/beam 
@ 1.2 mm) and the brightest HMPOs not detected with 
CARMA (- 250 mJy/beam). Thus the lack of CARMA 
detections around sources with low peak 1.2 mm flux den- 
sities, but otherwise large integrated 1.2 mm flux densi- 
ties and far- infrared luminosities, may indicate that a 
compact disk or envelope has dissipated around these 
sources. 

Since the detection rate of 3 mm sources is not strongly 
biased by the range of distances, masses and luminosi- 
ties in our sample, we investigate if the evolutionary 
stage correlates with the detection rate of compact 3 mm 
sources. The signposts of the age sequence that we are 
testing, in order of increasing age, are: (1) presence of 
6.7 GHz CH3OH and/or 22 GHz H2O maser emission and 
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no cm continuum emission, (2) both maser and cm con- 
tinuum emission, and (3) cm continuum emission but 
an absence of masers (iLo et al.111 975': 'G enzel fc Downesl 
[19771 iCodella et afl 119941 : iCode lla & Fellil 119951) . Fbr 
each of the 14 massive cores in our survey, we hst the 
presence of 3 mm continuum, cm continuum, and maser 
emission in Table [5] 

We detect 3 mm continuum emission from 3/5 of the 
"youngest" HMPOs, 4/5 of the "intermediate" aged HM- 
POs and 0/4 of the "oldest" HMPOs. Thus, the 3 mm 
continuum emission is preferentially detected toward the 
sources with maser emission, consistent with the pro- 
posed evolutionary sequence for high mass stars. Fur- 
ther, we would expect that the H2O masers and the 
6.7 GHz CH3OH masers to be nearly coincident with 
the compact 3 mm sources, since VLBI observations sug- 
gest that these masers originate from shocks associated 
with winds or jets near the base of the mol ecular outflow 
(jGoddi et al.l 120051 : iMoscadelh et al.ll2005D . All thirteen 
H2O masers with interferometric positions, and five of 
the eight 6.7 GHz CH3OH masers, are in fact located 
within 1.5" of a compact 3 mm continuum source. It is 
possible that the CH3OH maser sources without associ- 
ated 3 mm continuum emission may have too low of disk 
or envelope mass to be detected in our observations. 

In summary, all of the H2O and most of the 6.7 GHz 
CH3OH maser emission appear to be related to mas- 
sive disks or envelopes, but these disks/envelopes are 
not present toward sources that have radio continuum 
emission but no maser emission. The probability of as- 
sociation between HMPOs with 3 mm continuum emis- 
sion and HMPOs with 22 GHz H2O masers and 6.7 GHz 
CH3OH masers is 98% and 71%, respectively. These 
3 mm continuum sources are also likely to be driving 
outflows, and there is a 98% likelihood that 3 mm con- 
tinuum emission is associated with two possible outflow 
tracers, SiO (2-1) emission and 95 GHz CH3OH masers. 
There is a 99.8% likelihood of association between these 
two spectral lines. These observations are consistent with 
the hypothesis that as an H II region around a high mass 
protostar expands, it destroys the surrounding disk and 
envelope, and consequently shuts off the source of the 
outflows and H2O and CH3OH masers. 

6. SUMMARY 

In this paper we tested the hypothesis that H2O and 
CH3OH masers are signposts of an early phase in the 
evolution of a high mass protostellar object before an ex- 
panding UCHII region has destroyed the accretion disk. 
We used CARMA to observe the 3 mm continuum emis- 
sion at 1" resolution around a sample of HMPOs with 
and without masers and cm continuum emission. We also 
observed the 95 GHz CH3OH line to determine its asso- 
ciation with the thermal emission from dust surrounding 
the protostars. 

We detect (at > 3ct significance) 3 mm continuum emis- 



sion from all seven HMPOs with MAMBO 1.2 mm con- 
tinuum peak fluxes (taken from Beuther et al. 2002a) 
above 300 mJy/(ll")2, while none of the seven HMPOs 
with peak fluxes below this threshold were detected in 
our survey. We argue that this is not merely an observa- 
tional selection effect, because all 14 HMPO candidates 
could have been detected in the CARMA maps if their 
1.2 mm peak fluxes originated from disks or compact en- 
velopes, and were therefore concentrated into the ^1" 
CARMA beam. The detection rate of HMPOs at 3 mm 
is not correlated with distance to the HMPO, the 1.2 mm 
integrated flux detected with MAMBO, the total lumi- 
nosity of the HMPO or its mass (see Figure 2]). 

All seven HMPOs detected at 3 mm are associated with 
H2O masers, and only two of the seven HMPOs not de- 
tected at 3 mm are associated with H2O masers. The 
22 GHz H2O masers and 95 GHz CH3OH lines are signif- 
icantly (98% probability) associated with the detection 
rate of 3 mm continuum emission. There is no significant 
association between the detection rate of 3 mm contin- 
uum and 6.7 GHz CH3OH or 1.7 GHz OH masers or cm 
continuum emi ssion. The presence of S iO (2-1) emission 
in the study of iSridharan et ahl (|2002f ) was also a useful 
predictor of CARMA 3 mm continuum detections (98% 
probability) and 95 GHz CH3OH masers (99.8% proba- 
bility). 

The general picture of the evolution of high-mass 
protostars is consistent with our observations that the 
3 mm continuum emission comes from a disk and/or 
compact envelope around an HMPO. Our observations 
suggest that H2O masers are associated with these 
disks/envelopes, and that they are present both prior 
to and during the formation of an UCHII region. When 
the UCHII region disrupts the disk and compact enve- 
lope around the HMPO, the H2O and 95 GHz CH3OH 
masers are likely to be cut off and the 3 mm continuum 
emission falls below our threshold of detectability. 
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TABLE 1 

Source Properties 



TD A O 1\T ^ 1^-. „ 


1 

Group 


distance 


\r 2 
Vlsh 


log(Ltotal/L0) 


Feak Oiy(1.2mmJ 


Integratea 


mm) 


Core Masa 






(kpc) 


( km s ^ ) 




(mjy/beam) 


(Jy) 




(M©) 




I 


9.5 


95.9 


4.6 


265 




1.4 


6860 


1 8440—01 48 




8.3 


97.6 


4.7 


120 




0.5 


1717 


18517+0437 


1 


2.9 


43.9 


4.1 


812 




7.2 


2310 


18566+0408 


1 


6.7 


85.2 


4.8 


486 




1.6 


2110 


23151+5912 


1 


5.7 


-54.4 


5.0 


406 




2.0 


1229 


19217+1651 


2 


10.5 


3.5 


4.9 


640 




2.6 


9518 


20126+4104 


2 


1.7 


-3.8 


3.9 


1087 




5.8 


460 


20332+4124 


2 


3.9 


-2.0 


4.4 


265 




3.5 


1529 


23033+5951 


2 


3.5 


-53.1 


4.0 


631 




3.5 


2327 


23139+5939 


2 


4.8 


-44.7 


4.4 


530 




2.3 


1759 


19220+1432 


3 


5.5 


68.8 


4.4 


256 




1.9 


3406 


20205+3948 


3 


4.5 


-1.7 


4.5 


104 




0.9 


548 


22134+5834 


3 


2.6 


-18.3 


4.1 


229 




2.5 


436 


22570+5912 


3 


5.1 


-46.7 


4.7 


215 




2.1 


1469 



^ Group 1: With 22 GHz H2O and/or 6.7 CH3OH mascis and 110 cm continuum, Group 2: With 22 GHz H2O and/or 6.7 CH3O H mascrs 
and cm continuum, Group 3: Without 22 GHz H2O or 6.7 CH3OH mascrs and with cm continuum^ From [Sridharan" et al.l 1 120021) ^ 1.2 mm 
continuum flux densities fromlB^cuthcr ct all l|2002aD . measured in a ll" bcam'^ Mass of gas and dust inferred from 1.2 mm continuum flux 
densities, taken from lSeuther et al.| (1200211) 
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TABLE 2 

CARMA Observing Parameters 



IRAb Name 




RA 


1 

Dec 




Phase Calibrators 


Resolution 


Continuum Noise 


Noise/ channel'^ 






J2000 


J2000 






(" X ") 


(mjy/beam) 


(mJy/beam) 




18: 


37:17 02 


— 06-38-3n 


.70 


17tS1+n96 1743 — 038 


1.6x0.7 


1.3 


62 


1 8440—01 48 


18: 


■46-36 56 


—01:45:21 


.20 


17c;i_i_nQfi 1743 — 038 


1.6x0.7 


0.9 


34 


18517+0437 


18: 


:54:14.32 


+04:41:39, 


.69 


1751+096,1743-038 


1.4x0.7 


1.2 


49 


18566+0408 


18: 


:59:10.02 


+04:12:14, 


.69 


1751+096,1743-038 


1.6x0.7 


0.9 


26 


19217+1651 


19: 


:23:58.77 


+16:57:44, 


.80 


1751+096,1925+211 


1.5x0.7 


1.2 


31 


19220+1432 


19: 


:24:20.05 


+14:38:03, 


.60 


1751+096,1925+211 


1.4x0.7 


1.0 


32 


20126+4104 


20: 


:14:25.86 


+41:13:33, 


.99 


BLLAC, 2007+404 


1.4x0.8 


1.0 


37 


20205+3948 


20: 


:22:20.87 


+39:58:15, 


.00 


BLLAC, 2007+404 


1.6x0.8 


1.5 


78 


20332+4124 


20: 


:34:59.72 


+41:34:49, 


.40 


BLLAC, 2007+404 


1.2x0.6 


1.0 


39 


22134+5834 


22: 


:15:09.51 


+58:49:05, 


.99 


BLLAC, 2007+404 


1.5x0.8 


1.3 


61 


22570+5912 


22: 


:59:05.37 


+59:28:19, 


.18 


BLLAC, 0102+584 


1.0x0.8 


0.6 


30 


23033+5951 


23: 


:05:25.31 


+60:08:06, 


.28 


BLLAC, 0102+584 


1.0x0.8 


0.6 


22 


23139+5939 


23: 


:16:10.45 


+59:55:28, 


.48 


BLLAC, 0102+584 


1.0x0.8 


0.6 


26 


23151+5912 


23: 


:17:21.02 


+59:28:48, 


.48 


BLLAC, 0102+584 


1.2x0.7 


0.8 


34 



Phase center positions taken from 1.! 
natural weig hting^ Channel width of 0.488 MHz 



TABLE 3 

Properties of Detected 3 mm Continuum Sources 



IRAS Name 


RA 


Dee 


(Aa, A5)i 






/d 2 


Mass 


Density 




J2000 


J2000 


("/') 


(mJy/beam) 


(mJy) 




(Mo) 


(10^ cm-3 ) 


18517+0437 


18:54:14.24 


+04:41:40.7 


(-1.1, 1.0) 


8.5+1.4 


16+2 


1.0 


74 


8.6 


18566+0408 


18:59:09.99 


+04:12:15.3 


(-0.4, 0.6) 


4.4+1.2 


10+2 


1.0 


250 


2.3 


19217+1651 


19:23:58.81 


+16:57:41.1 


( 0.6,-3.7) 


55+3 


66+1 


0.3 


1200 


2.9 


19217+1651 


19:23:58.70 


+16:57:42.0 


(-1.1,-2.8) 


4.9+1.3 


9+2 


1.0 


540 


1.3 


20126+4104 


20:14:26.03 


+41:13:32.6 


( 1.8,-1.4) 


18±1 


33+2 


0.9 


47 


27 


23033+5951 


23:05:24.65 


+60:08:09.3 


(-4.9, 3.0) 


4.4+0.9 


8+1 


1.0 


54 


3.6 


23033+5951 


23:05:24.97 


+60:08:14.2 


(-2.5, 7.9) 


2.3+0.6 


6+1 


1.0 


40 


2.7 


23033+5951 


23:05:25.09 


+60:08:16.3 


(-1.7,10.0) 


3.4+0.9 


8+1 


0.4 


21 


1.4 


23139+5939 


23:16:10.45 


+59:55:28.6 


( 0.0, 0.1) 


2.8+0.5 


12+2 


0.7 


110 


2.7 


23151+5912 


23:17:20.94 


+59:28:47.6 


(-0.8,-0.9) 


3.5+1.2 


8+2 


1.0 


140 


2.2 



Difference between the measured 3 mm continuum source and the 1.2 mm core (see Table [TJ^ Fraction of 3 mm 
emission from dust (see Section [4t 



TABLE 4 

Properties of Detected 95 GHz CH3OH Spectral Lines 



IRAS Name 


RA 


Dec 




separation^ 


'5'pcak 


Vlsr 


FWHM 






J2000 


J2000 




(10-* AU) 


(Jy) 


(kms^^) 


(kms~^) 


18517+0437 


(a) 


18:54:14.74 


+04:41: 


:42.6 


2.22+0.05 


2.98 


43.4 


1.2 


18517+0437 


(b) 


18:54:14.46 


+04:41: 


:44.6 


1.48+0.05 


1.60 


43.9 


1.4 


20126+4104 


(a) 


20:14:26.72 


+41:13: 


:29.8 


1.40+0.01 


2.72 


-4.2 


1.5 


20126+4104 


(b) 


20:14:25.24 


+41:13: 


:34.9 


1.57+0.01 


1.46 


-2.1 


2.2 


20126+4104 


(c) 


20:14:25.16 


+41:13: 


:36.4 


1.79+0.01 


2.92 


-2.4 


2.0 


20126+4104 


(d) 


20:14:25.44 


+41:13: 


:37.3 


1.39+0.01 


6.14 


-2.9 


1.4 


20126+4104 


(e) 


20:14:25.41 


+41:13: 


37.9 


1.49+0.01 


6.35 


-2.9 


1.3 


23033+5951 


23:05:24.59 


+60:08: 


:09.4 


0.16+0.04 


4.55 


-54.3 


1.4 


23139+5939 


(a) 


23:16:10.86 


+59:55: 


:20.7 


4.07+0.13 


0.44 


-44.6 


1.7 


23139+5939 


(b) 


23:16:10.40 


+59:55: 


:28.0 


0.33+0.13 


1.23 


-45.9 


1.8 


23151+5912 


23:17:21.90 


+59:28: 


:45.7 


4.30+0.22 


1.87 


-52.8 


1.4 



Projected separation between 95 GHz CH3OH maser position and nearest 3 mm continuum peak. 
Uncertainties estimated from signal to noise of the maser and continuum detections. 
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TABLE 5 

Association of Masers and Continuum Detections 



IRAS Name 


Group 


3 mm Continuum 


3.6 cm Continuum'^ 


H2O Maser'^ 


CH3OH Maser'^ 


CH3OH Maser^ 


OH Maser* 










22 GHz 


6.7 GHz 


95 GHz 


1.7 GHz 


1 R'^d'^—(^(^A^ 

iOOIiJ — UU^ -L 




N 


N 


Y 


Y 


N 


Y 






N 


N 


N 

ly 


Y 


N 

ly 


Y 


18517+0437 


1 


Y 


N 


Y 


Y 


Y 


N 


18566+0408 


1 


Y 


N 


Y 


Y 


N 


Y 


23151+5912 


1 


Y 


N 


Y 


N 


Y 


N 


19217+1651 


2 


Y 


Y 


Y 


Y 


N 


Y 


20126+4104 


2 


Y 


Y 


Y 


Y 


Y 


Y 


20332+4124 


2 


N 


Y 


Y 


N 


N 


N 


23033+5951 


2 


Y 


Y 


Y 


N 


Y 


N 


23139+5939 


2 


Y 


Y 


Y 


Y 


Y 


Y 


19220+1432 


3 


N 


Y 


N 


N 


N 


Y 


20205+3948 


3 


N 


Y 


N 


N 


N 


N 


22134+5834 


3 


N 


Y 


N 


N 


N 


N 


22570+5912 


3 


N 


Y 


N 


N 


N 


N 



^ Group 1: With 22 GHz H2O and/or 6.7 CH3OH masers and no cm continuum. Group 2: With 22 GHz H2O and/or 6.7 CH3OH masers and 
cm continuum, Group 3: Without 22 GHz H2O or 6.7 CH3OH masers and with cm continuum^ This studv^ iSridliaran et"al] | |2002D "^ lEdris et al.l 
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Fig. 1. — Contour plots of the 3 mm continuum flux density for HMPO candidates sources detected with CARMA at the > 3cr level, 
centered at the position of the peak 1.2 mm MAMBO (11" resolution) emission. Grey contours show emission at 3, 6 and 12 times the 
noise in each map, and black contours show emission at —3, —6 and —12 times the noise. Red crosses and triangles show the positions of 
22 GHz H2O and 6.7 G Hz CH3OH masers, respectively, and blue 'X' mark the positions of 3.6 cm continuum sources (Bcuthcr ct al. 200251; 
ISridharan et al1l2002l . private communication with Beuther). The source IRAS 18517-1-0437 was not observed by Bcuthcr ct al. (2002b|). 
Green squares show the 95 GHz CH3OH maser emission detected by CARMA. The resolution of the 3 mm continuum images is 
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Fig. 2. — Contour plots of the 3 mm continuum flux density for eight HMPO candidates not detected with CARMA. Grey contours 
show fluxes at 2.0 and 2.5 times the noise in each map, and black contours show fluxes —2.5 and —2.0 times the noise. The symbols are 
otherwise the same as presented in Fig. [l] 
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Fig. 3. — Spectra of the 95 GHz CH3OH maser detections, with the same identifications shown in Figure[T]and Table|4] 
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Fig. 4. — Properties of the HMPO candidates with (black circles) and without (gray circles) CARMA 3 mm continuum detections. 
Top: RMS noise in the CARMA 3 m m im age versus the distance to the HMPO. Middle: Mass of the dense core, as measured by 1.2 mm 
continuum images ( Beuther et al.l;2002al). versus the bolometric luminosity. Bottom: Integrated flux density in the 1.2 mm single dish 
maps versus the peak 1.2 mm flux density. These figures show that HMPOs with and without CARMA detections have similar noise 
characteristics, distances, core masses, and luminosities. However, sources without CARMA continuum detections tend to have lower peak 
flux densities. The data displayed in these plots are presented in Tables [T] and [2] 



